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The PowerFLOW® code has been used to perform numerical simulations of the high-lift version of the Com- 
mon Research Model (HL-CRM) that will be used for experimental testing of airframe noise. Time-averaged sur- 
face pressure results from PowerFLOW® are found to be in reasonable agreement with those from steady-state 
computations using FUN3D. Surface pressure fluctuations are highest around the slat break and nacelle/pylon 
region, and synthetic array beamforming results also indicate that this region is the dominant noise source on 
the model. The gap between the slat and pylon on the HL-CRM is not realistic for modern aircraft, and most 
nacelles include a chine that is absent in the baseline model. To account for those effects, additional simulations 
were completed with a chine and with the slat extended into the pylon. The case with the chine was nearly 
identical to the baseline, and the slat extension resulted in higher surface pressure fluctuations but slightly re- 
duced radiated noise. The full-span slat geometry without the nacelle/pylon was also simulated and found to be 
around 10 dB quieter than the baseline over almost the entire frequency range. The current simulations are still 
considered preliminary as changes in the radiated acoustics are still being observed with grid refinement, and 
additional simulations with finer grids are planned. 


Nomenclature 
Cp coefficient of pressure Supenscrin : ; ; 
é local deployed chord perturbation quantity (e.g., p = p — po) 
pressure Subscript: 
rms root mean square & dimensional reference quantity 
VR variable resolution 


L,Y, z Cartesian coordinates 


I. Introduction 


Aircraft noise reduction, including that of the airframe, is an important goal of the NASA Advanced Air Transport 
Technology (AATT) Project, which is supporting a combined experimental and computational effort to better understand 
and mitigate the sources associated with slat noise. The nonpropulsive (or airframe) sources of aircraft noise include 
high-lift devices (e.g., the leading-edge slat and trailing-edge flaps) and the aircraft undercarriage. The ranking of these 
sources is configuration dependent; however, both model-scale tests!’ and flyover noise measurements® have identified 
the leading-edge slat as a prominent source of airframe noise during aircraft approach. To further develop airframe 
noise reduction technology, NASA is currently planning to construct a 10%-scale version of the High-Lift Common 
Research Model (HL-CRM) recently developed by Lacy and Sclafani.? The original cruise configuration NASA CRM is 
an open geometry that has been widely used in the AIAA Drag Prediction Workshops.'° The NASA CRM! consists of a 
contemporary supercritical transonic wing with flow-through nacelles and a fuselage that is representative of a widebody 
commercial transport aircraft. The new HL-CRM is also an open geometry that is being used in the AIAA Geometry 
and Mesh Generation Workshop!” and 3 AIAA High-Lift Prediction Workshop.!* 
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Two views of the HL-CRM are shown in Fig. 1. The geometry includes inboard and outboard flaps that meet in 
the center. There are also inboard and outboard slats, but there is a gap between them to accommodate the pylon for a 
flow-through nacelle. In the landing configuration, both flap deflections are set to 37° and the slats are set to 30°. The 
current geometry does not have any brackets nor some of the modifications that will be necessary during the detailed 
design of the wind tunnel model. However, simulations of the basic geometry should still give an initial indication of 
both aerodynamic and acoustic performance of the HL-CRM. 

NASA is currently developing an Active Flow Control version of the HL-CRM, and both the conventional and flow 
control semispan models will be tested in the NASA Langley Research Center (LaRC) 14- by 22-foot (14x22) subsonic 
tunnel in 2018/2019. The entry will include aeroacoustic measurements using a microphone array!* that was previously 
used during the testing of an 18% Gulfstream Aircraft. During that campaign, various flap and landing gear noise 
reduction devices!> were developed and evaluated. However, that model did not include a slat, so no slat noise reduction 
devices were developed. To demonstrate an overall reduction in airframe noise for large commercial transport aircraft, 
the HL-CRM will be used as a platform to evaluate slat noise-reduction concepts such as the slat-cove filler'® '7 and 
slat-gap filler'® at a higher technology readiness level. Slat-cove fillers were tested on a trapezoidal-wing model!® and 
the 26% 777 STAR model, but those treatments were incapable of being stowed. The slat gap- and cove-fillers that will 
be tested on the HL-CRM will be constructed out of shape-memory alloys so that the slat could be stowed. However, the 
HL-CRM slat will not articulate, and other testing will be used to evaluate additional structural aspects of the designs. 

Computational simulations are being used to support the model development and to aid in the design of the noise 
reduction devices. Although several computational fluid dynamics (CFD) codes are being employed in the overall effort, 
the commercial CFD software PowerFLOW® version 5.3c is being used to make aeroacoustic predictions of the noise 
from the HL-CRM. PowerFLOW® was used extensively during the design of the noise reduction technology applied 
the to Gulfstream aircraft model tested in the LaRC 14x22,”°?! and the noise predictions made before the experimental 
testing compared very well with the measurements. PowerFLOW® was also used for slat noise simulations involving the 
unswept 30P30N high-lift configuration from the BANC series of workshops.”* Preliminary simulations of the HL-CRM 
in the landing configuration have been completed with PowerFLOW® and the mean flow field will be shown to be in 
reasonable agreement with the steady CFD results” from the FUN3D code.”* Surface pressure fluctuations and synthetic 
microphone array beamform maps will be used to identify potential noise sources. Some unexpectedly prominent noise 
sources have been identified in the vicinity of the nacelle/pylon junction where there is a break between the inboard and 
outboard slat sections. Several means of mitigating these junction related noise sources are then evaluated numerically, 
but only the complete elimination of the nacelle/pylon was found to allow the normal slat and flap sources to be easily 
identified. 


II. Simulation Methodology 


The numerical simulations presented in this paper were performed using the commercial CFD software PowerFLOW® 
which is a compressible flow solver based originally on the three-dimensional 19 state (D3Q19) Lattice Boltzmann 
Model (LBM). The PowerFLOW® code represents LBM-based CFD technology developed over the last 30 years,>-° 
and has been extensively validated for a wide variety of applications ranging from academic direct numerical simula- 
tions (DNS) cases to industrial flow problems in the fields of aerodynamics* and aeroacoustics.*! In contrast to methods 
based on the Navier-Stokes (N-S) equations, LBM uses a simpler and more general physics formulation at the micro- 
scopic level.> The LBM equations recover the macroscopic hydrodynamics of the Navier-Stokes equations*~ ** through 
the Chapman-Enskog expansion. The local formulation of the LBM equations allows a highly efficient implementation 
for distributed computations on thousands of processors. The low dissipation and dispersion properties of the numerical 
scheme produces aerodynamic and aeroacoustic results that are generally comparable to large eddy simulations obtained 
with classical CFD solvers, as shown in Refs. 34 and 35, and demonstrated in the comparative study of flow over tandem 
cylinders by Lockard.*° 

The classical LBM scheme is typically valid in the low speed regime up to a local Mach number of 0.5. Recent 
extensions of the scheme recover a fully unsteady compressible form of the Navier-Stokes equations.*’~*° Applications 
of this new version at transonic conditions were presented in recent papers by Koeing and Fares‘? for the NASA CRM, 
and by Duda et al.*! for a sweeping jet (fluidic) actuator operating at choked conditions. The newer version of the 
PowerFLOW® code, which incorporates the modified LBM scheme suitable for simulating flows containing transonic 
flow regimes, has two options for higher-speed flows: the high-subsonic option for 0.5 < Mach < 0.9, and a transonic 
option for 0.9 < Mach < 2.0. All of the results presented here were obtained using the high-subsonic option of the 
baseline solver. The local maximum time-averaged Mach number for the landing conditions under investigation is 0.6, 
which is low enough for the baseline solver. However, some transients will exceed this value, so the transonic version 
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(not used here) may produce somewhat improved results. 

The PowerFLOW® code can be used to solve the Lattice-Boltzmann equation in a DNS mode,*” where all of 
the turbulent scales are spatially and temporally resolved. However, for most engineering problems at high Reynolds 
numbers, the simulations are usually conducted in conjunction with a hybrid turbulence modeling approach where the 
small scales are modeled and the largest scales containing most of the energy are directly resolved. The current work the 
Lattice-Boltzmann Very Large Eddy Simulation (LB-VLES) approach described in Refs. 27, 28 and 43 is used. 

The standard Lattice-Boltzmann boundary condition for no-slip or the specular reflection for free slip condition are 
generalized through a volumetric formulation?>: *° near the wall for arbitrarily oriented surface elements (surfels) within 
the Cartesian volume elements (voxels). This formulation of the boundary condition on a curved surface cutting the 
Cartesian grid automatically conserves mass, momentum, and energy, and is compatible with the general second-order 
spatial accuracy of the underlying LBM numerical scheme. To reduce the resolution requirement near solid surfaces for 
high Reynolds number flows, a hybrid wall function is used to model the near wall region of the boundary layer.*? “4 

The Lattice-Boltzmann equation is solved on embedded Cartesian meshes, which are generated automatically within 
the flow solver on the basis of input specifications provided by the user. Variable resolution (VR) regions can be defined 
to allow for local mesh refinement of the grid by successive factors of two in each direction.” The PowerFLOW® code 
scales well on modern computer clusters consisting of thousands of processors, making it ideally suitable for large scale 
applications. The LBM methodology described here has been extensively validated for a wide variety of applications 
ranging from DNS for academic cases”? to LB-VLES for industrial flow problems in the fields of aerodynamics, thermal 
management, and aeroacoustics, see e.g., Refs. 30 and 44-47. 


III. Results 


PowerFLOW® simulations of the semispan HL-CRM have been completed on a series of meshes. A planar view of 
a typical mesh is shown in Fig. 2. The increase in resolution near solid surfaces is obvious, but finer resolution was also 
specified to capture important flow features in the interior of the flow-field. Figure 3 shows some of the VR regions on 
the top and bottom surfaces of the wing. The wireframe objects in pink have the finest resolution. Some of these regions 
are boxes and cylinders defined to enclose important flow features. In addition, the region in the slat cove is encompassed 
by a surface that was defined based on an isocontour obtained from a preliminary solution. The grid spacing one level 
down, which represents a coarsening by a factor of two in each coordinate direction, is shown by the red outlined objects. 
An additional level down is shown by the green outlined objects. These regions were defined by examining solutions 
on coarser grids. Important features were identified from contours of the steady and unsteady surface pressure, surface 
and volume streamlines, and isosurfaces. Furthermore, noise source regions identified by synthetic microphone array 
analyses were also targeted for refinement. 

The first set of grids was denoted as version zero (vO), with both coarse (c) and medium (m) spacings for the finest 
cells. The smallest grid, vOc, had approximately 208 million voxels and a minimum grid spacing of 0.432 mm. Most of 
the refinement at level 0 was dictated by the distance from solid surfaces. The vOm grid was uniformly refined by a factor 
of 1.5. The version | grids primarily included refinement based on location such as around leading and trailing edges as 
well as global refinement. The vic grid had a minimum spacing of 0.216 mm, and 0.144 mm for vlm. The version 2 
grids had more targeted refinement around flow features, and the v2c grid had 900 million voxels and the same minimum 
grid spacing of 0.216 mm as the vic grid. The meshes are summarized in Table | where the number of Fine Equivalent 
(FE) Voxels is also included. PowerFLOW® only updates cells in time as needed based on the size of the cell, and the 
Fine Equivalent is an estimate of the average number of cells that must be updated at each time step. However, when 
comparing grids with different minimum spacings, the number of time steps required to reach a fixed point in time will 
vary linearly by the size of the smallest voxel. Localized refinement around the pylon and wing tip is shown in Fig. 4. 
The version 2 grids have larger regions of fine grid in locations where strong vorticity and unsteadiness is likely to occur. 
However, even much finer grids may be necessary to capture all of the relevant noise generating mechanisms. 

All of the simulations have been run at landing conditions with a Mach number of 0.20, static temperature of 15°C 
(519°R), Reynolds number (based on the mean aerodynamic chord of 0.7 m or 27.58 in) of 3.27x10® and an angle 
of attack of 8° This Reynolds number corresponds to the conditions expected for the 10% model that will be tested in 
the NASA LaRC 14x22 tunnel. The wing semispan is 2.938 m (115.675 in), which corresponds to a large transport 
aircraft. The slats are deployed at an angle of 30° and the flaps at 37° The mean aerodynamic quantities such as lift 
and drag did not vary significantly across any of the meshes used in PowerFLOW® and they compare reasonably well 
with steady-state results’? from FUN3D. The FUN3D Reynolds-Averaged-Navier-Stokes calculations were obtained 
with the Spalart-Allmaras*® turbulence model on a grid with 92 million nodes. No grid refinement study was performed 
with FUN3D, but the grid is considered relatively fine for this type of configuration. The FUN3D solutions specified a 
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Table 1. Summary of Grid Parameters. 


ID Mindx(mm)_ Voxels (x10°) FE Voxels (x 10°) 


vOc 0.432 208 47 

vOm 0.288 500 100 
vic 0.216 362 89.5 
vim 0.144 916 211 
v2c 0.216 900 338 


symmetry condition at the centerline of the model in an otherwise free-air setup, whereas the PowerFLOW® simulations 
had the semispan model, with an 8.89 cm (3.5 in) extension that would be required for wind-tunnel testing, mounted 
on a viscous flat plate. Hence, there are some differences in the configurations. PowerFLOW® has also been used to 
perform simulations of the model in the LaRC 14x22 wind tunnel operating in its closed mode with walls enclosing the 
test section. The lift and drag from an angle of attack sweep was compared with the results from calculations with the 
tunnel in the open mode with the top and side walls removed. The lift curve slopes were similar, but the closed wall 
all the had higher lift at the same angle of attack. For the landing configuration under investigation, an angle of attack 
adjustment is sufficient to obtain an aerodynamic equivalence. Therefore, because reflections from tunnel walls present 
difficulties for acoustic calculations, the free-air or open tunnel mode is used for all the simulations. 


IV. Time-Averaged Pressure Results 


A comparison of time-averaged coefficient of pressure, C;,, contours is shown in Fig. 5. On the top or suction side, 
the primary difference is around the outboard portion of the flap where FUN3D predicts separation, but the flow remains 
mostly attached with PowerFLOW® on this v2m grid as well as the coarser grids. The comparison is similarly good on 
the bottom or pressure side, with most of the differences being around the pylon or its wake. 

Close-up comparisons of C,, are shown in Fig. 6, with the pylon region presented in (a) and (c), where the agreement 
is reasonably good. However, the footprint of the flap-side-edge vortex is clearly visible in the PowerFLOW® result 
in Fig. 6 (d), but it is absent in the FUN3D solution in (b). Although some discrepancies are expected because of 
the different numerical approaches in the methods, the discrepancy around the flap is likely related to the separation 
on the outboard flap predicted by FUN3D. Until the experiments are completed, the correct state of the flow over the 
flap suction surface will be uncertain. However, additional information may be available from the 3"! AIAA High-Lift 
Prediction Workshop!* where many codes will be used to simulate the HL-CRM. Unfortunately, the configuration being 
investigated for the workshop does not have the nacelle/pylon installed. 

More quantitative comparisons of the surface pressure are made along the planar cuts in Fig. 7. The cuts are indicated 
by the colored lines on the wing, and the parameter 7 is the distance in the spanwise direction normalized by the total 
semispan length. The comparisons are shown in Fig. 8 with the most inboard cut at 7 = 0.240 in (a) and the most 
outboard cut at 7 = 0.997 in (f). The streamwise coordinate has been normalized by the local deployed chord length. 
The PowerFLOW® results from six different meshes are included, with the mesh parameters given in Table 1. In general, 
these time-averaged PowerFLOW® results are insensitive to the grid resolution. The lone exception is near the wing tip 
where there is some small variation in the wake of slat tip in Fig. 8 (e) and near the trailing edge of the wing tip (f). The 
agreement between FUN3D and PowerFLOW® is quite good at the two inboard stations, but degrades on the outboard 
flap region where FUN3D predictions indicate flap separation. At the two most outboard stations, the agreement is again 
reasonable except in the same regions where some variation is seen with mesh resolution in the PowerFLOW® results. 


V. Fluctuating Pressure Results 


FUN3D has only been run in a steady-state mode, so no unsteady data is available for comparison. Therefore, only 
PowerFLOW® results will be presented in the remainder of the paper. The contours of the root-mean-square (rms) of the 
fluctuating surface pressure from PowerFLOW® are shown in Fig. 9. In both the top and bottom views, the fluctuations 
around the nacelle/pylon are the highest and most extensive. The unexpectedly high surface pressure fluctuation levels 
around the nacelle prompted an examination of the pylon/wing juncture region on modern aircraft, which revealed that 
the gap between the inboard slat edge and the pylon is usually filled by some device. Images of a Boeing 787 and an 
Airbus A320 in the landing configuration are shown in Fig. 10. A comparison of the pylon region in Fig. 11 reveals that 
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an important difference between the configurations involves the presence of a Krueger flap between the inboard slat and 
the pylon on the Boeing 787 and a bulge in the pylon on the A320 that closes the gap with the inboard slat. Curiously, 
we have not found any aircraft where the gap between the outboard slat and the pylon is closed in a similar fashion. 
Another obvious difference is the relatively narrow and more rectangular pylon on the HL-CRM in comparison with 
the two real aircraft. The pylon and gap treatments on real aircraft are most likely designed for aerodynamic reasons, 
but they are also likely to have an important effect on the acoustics. One of the goals of the upcoming HL-CRM test is 
to examine slat-cove noise reduction devices, and minimizing other sources, especially those that are unrealistic, which 
will make it easier to isolate true slat-cove noise. Therefore, additional PowerFLOW® calculations have been performed 
with a nacelle chine, with the slat extended into the pylon, and with a full-span slat and no nacelle/pylon. All of these 
simulations were run with meshes similar to the v2c grid for the baseline geometry. The nacelle chine was taken from 
another aircraft geometry of a similar sized aircraft, but may not be identical to the chine that is ultimately designed for 
the HL-CRM. Furthermore, only three locations of the chine were evaluated at an angle of attack of 8? and the position 
where the vortex coming off the chine most closely passed through the slat break region was chosen for an additional 
simulation. However, this may not be an optimal position for the stall case, i.e., the primary target for the employment 
of a chine. 

The rms of the perturbation pressure from the four simulations is shown in Fig. 12. The nacelle and pylon have 
been removed to allow a direct view of the slat break region. The chine has very little effect on the fluctuations, and 
the extended slat has more pronounced fluctuations on the inboard slat near the slat break. As should be expected, the 
full-span slat case has much lower fluctuation levels on the slat. All of the simulations exhibit significant unsteadiness 
near the wing-fuselage junction, but the current version of the HL-CRM does not have a slat horn or fairing that would 
have smoothed out the flow in this region. Fluctuations are also seen around the inboard flap edge and the fuselage. 
Some of this is caused by the flow around the flap edge, but some is because of the fuselage fairing, which is not readily 
apparent in the figure. The fuselage fairing is not very well resolved and some of the unsteadiness may disappear with 
grid refinement in this region. Additional regions with significant fluctuation levels include the outboard flap tip and the 
main wing cove, especially in the wake of the pylon. The fluctuation levels in the main cove were observed to increase 
significantly after refining the grid in the wake of the slat edges adjacent to the pylon. 

Close-up views of the nacelle/pylon are presented in Fig. 13 where now the case with the chine exhibits some subtle 
differences with the baseline. Surface streamlines have been superimposed on the contours to give another indication of 
the flow patterns. Again, the extended slat case appears to have enhanced fluctuation levels both on the slat and nacelle. 
The extended slat presents more blockage to the flow, so the fluid has to accelerate more to get around the pylon. The 
higher velocities may be partially responsible for the higher fluctuations. Another view of this region is presented in 
Fig. 14 where the complicated flow patterns and stagnation regions around the pylon are readily apparent in the surface 
streamlines. In contrast, the flow around the full-span slat is fairly uniform in the span and does not exhibit any regions 
with higher fluctuation levels as those seen in the three other cases. 

The results from all the simulations were relatively similar around most of the model, so only the results from the 
baseline case are presented for some other important regions in Fig. 15. The top of the wing tip is shown in Fig. 15 
(a), where the wake of the slat tip is evident. Only a small portion of the fluctuations are really from the slat edge as 
considerable separation also occurs on the edge created by the unslatted portion of the main element. After sufficient 
grid refinement, this vortical feature persisted so that its footprint was visible to the trailing edge of the wing. Significant 
fluctuations are also visible at the trailing edge near the wing tip. The wing tip vortex migrates from the side to the top 
at around the 80% chord position and is still relatively close to the wing surface as it moves across the trailing edge. The 
fluctuations near the leading edge of the wing tip are more visible in Fig. 15 (b) and are associated with disturbances in 
the slat cove being carried out past the slat edge by a strong spanwise flow and then being wrapped around the wing tip. 
An additional region with significant fluctuation levels includes the outboard flap edge shown in Fig. 15 (c). A typical 
dual vortex system is observed where the side-edge vortex migrates to the top around the 60% chord position and merges 
with the upper surface vortex. The last region examined is around the wing/fuselage junction shown in Fig. 9 (d) where 
high fluctuation levels are observed just upstream of the wing and around the edge created by the unslatted portion of 
the main element. 


VI. Acoustics: Synthetic Array Beamforming 


Although the simulations indicate strong unsteadiness in several locations on the model, the rms levels include 
fluctuations at all frequencies. Some very low frequency oscillations may be present in the pressure signals and may 
even dominate the rms levels. These are not particularly relevant for acoustics and, hence, contours focused on just the 
rms pressure fluctuations can be misleading. Even without this contamination, only a small fraction of fluctuating energy 
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is actually converted into acoustics. Hence, for any frequency, the rms pressure levels can only give an indication of 
potential noise source locations. One method to assess actual noise sources is through array beamforming. Conventional 
beamforming involves assuming an acoustic source basis function (such as a monopole in a uniform flow) and placing 
one of these sources at every point in a mesh surrounding a region where sound sources are expected. A source strength 
amplitude for each grid point is determined by how well cross-correlations of the signals across the microphone array 
are consistent with those based on the assumed basis function. Using the distances between the grid point and the 
microphones, each signal is adjusted in amplitude and time (or phase). They can then be combined such that the 
portion of the signal consistent with the assumed basis adds up constructively, whereas the inconsistent portion combines 
destructively. However, this assumes that the sources are uncorrelated and that their directivity is consistent with that 
of the basis function, which is typically a monopole. Furthermore, even when the source coincides exactly with the 
assumed basis, the array response is dependent on the particular arrangement of the microphones relative to the sources. 
Deconvolution methods have been developed to account for the array response that can provide spectra equivalent to what 
would be obtained by a single microphone, but all of these algorithms require certain assumptions and some of them can 
be computationally expensive. Nonetheless, microphone arrays have provided valuable information about noise sources 
when the elevated background noise would render single microphone measurements useless. In particular, the contour 
maps of source strength provide information about the location of sources that was not available previously. Array 
beamforming is typically used with experimental data, but the technique is now being applied to numerical simulation 
data? 

A methodology commonly used to make aeroacoustic predictions using CFD involves coupling the near-field solution 
from the CFD to an acoustic analogy such as the Ffowcs Williams and Hawking’s equation?! (FW-H). These predictions 
are often computed at the center of a microphone array and compared with the array output. However, with a minimal 
increase in computational cost, the predictions can be made at all microphone locations in an array, and the signals 
processed in the same manner as experimental data. Hence, the simulations provide synthetic array data that can be used 
with beamforming techniques. 

The far-field noise from the HL-CRM was calculated using the FW-H equation>!:** solver described by Bres.>* 
Results were obtained using 0.4 seconds of pressure time history data on all solid surfaces on the model but not on the 
flat plate. The array is positioned in the 90° position, geometrically directly beneath the center of the wing. The 97- 
microphone array is 276 inches (7 m) from the model centerline, with an outer diameter (microphone to microphone) of 
78.6 inches (2.0 m). An array shading algorithm was employed to exclude certain microphones based on the frequency 
so that sources appear similar in size across the frequency range and to reduce the distances between the included 
microphones as the frequency increases, which helps to minimize the detrimental effects of decorrelation. The array 
data was processed using the Exa beamforming code that uses the CLEAN-SC™ deconvolution approach. Detailed 
information about the array and processing procedure are in Refs. 14 and 50. 

CLEAN-SC beamform maps are shown in Fig. 16 for the four configurations under investigation. The data was 
processed in narrow band with a bin width of 125 Hz. The contour levels are normalized by the maximum for that 
map, and the range is 20 dB. Hence, only the relative distribution of source levels should be compared between the 
four cases. However, the baseline and chine cases are nearly identical with the extended slat being a few dB quieter for 
frequencies below | kHz and above 4 kHz. The full-span slat case is around 10 dB quieter than the baseline over almost 
all frequencies, with the exceptions being at 2 and 2.5 kHz where tones are present in the full-span slat spectra. As seen 
in Fig. 16, the source of the tone is associated with the outer section of the outboard slat. A similar source location was 
observed in the experimental testing of a 26% 777 model!®:>> and in a flight test.°° In the current simulations, the tones 
appear to be reminiscent of the narrow band peaks observed in many small-scale slat experiments and accompanying 
computations. The particular geometric arrangement and flow in the outboard section of the slat may lend itself to the 
narrow band peak phenonema. However, the current model does not have any brackets, which may alter the flow and 
have an impact on these tones. The three other cases in Fig. 16 are dominated by a source around the pylon, with the 
other sources being 10 dB down. Although the secondary sources are evident in these results, they may be difficult to 
observe in an experiment because of much higher background noise levels. The secondary noise sources include the 
wing tip, outboard flap edge, and slat break wake on the upper surface of the wing interacting with the main flap cove. 
This last source is the one seen on the inboard flap. The tonal slat source is not observed in the cases with the pylon 
because they do not have as much resolution in the slat cove. An additional simulation of the baseline with increased 
slat cove resolution has been completed, and the same slat tone phenomenon was evident. 

At a frequency of 5.0 kHz, Fig. 17 shows that the outboard flap edge source is now stronger than the one near the 
wing tip, but the pylon is still dominant. A source at the wing/fuselage junction is now more apparent as well. The maps 
for 20.0 kHz in Fig. 18 are indicative of those for the remainder of the frequency range up to 30 kHz. The pylon region 
persists as a strong source, but the wing tip and outboard flap are visible for the two quieter cases. Extensive sources are 
now visible on the fuselage for the full-span slat case, and we believe that these may be a result of insufficient resolution 
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around the fuselage fairing. 


VII. Conclusions 


An initial evaluation of the landing configuration of the HL-CRM has been completed using the PowerFLOW® 
code. Steady-state FUN3D results for the surface pressure compare reasonably well with time-averaged pressure from 
PowerFLOW® although FUN3D predicts much more extensive separation on the outboard flap. The effect of the grid 
density on the time-averaged quantities from the PowerFLOW® solutions was relatively small except in a few regions 
with strong vorticity and wakes such as at the wing tip. Although these regions only make up a small contribution to the 
overall lift and drag, they may have greater influence on the radiated noise. 

Surface pressure fluctuations were strongest around the slat break, nacelle/pylon region. A comparison of the HL- 
CRM geometry in this region with modern aircraft showed that the HL-CRM has a gap between the inboard slat edge 
and pylon that is usually closed by a Krueger flap or an enlarged pylon. Furthermore, a chine is usually installed on 
the nacelle. Simulations were run with a chine on the nacelle, the slat extended into the pylon, and with a full-span slat 
without the nacelle/pylon. Neither the chine nor the extended slat reduced the surface pressure fluctuations around the 
slat break. 

Noise source regions were identified using synthetic microphone analyses with deconvolution processing. The chine 
case was nearly identical to the baseline, and the extended slat resulted in a few dB of reduction over a portion of 
the frequency range of interest. A source around the wing/fuselage juncture was visible in all the cases, and further 
calculations with a slat horn are needed to see if such a fairing will reduce its intensity. Other sources were visible at 
the locations typically associate with airframe noise, such as the flap edges, wing tip, and slat cove. However, they were 
typically more than 10 dB weaker than the source at the pylon, and would be difficult to identify from experimental data 
with significant background noise. Only the full-span slat case was satisfactory for identifying all the typical airframe 
noise sources. 

Although several mesh densities have been used in the course of the simulations, the noise results are still varying 
too much to be confident in anything other than the trends. We are performing semispan simulations of a large transport 
aircraft at a relatively high Reynolds number with around | billion voxels (or 340 million fine equivalent voxels). To 
put things in perspective, PowerFLOW® simulations” attempting to resolve the slat noise sources from a 2-D extruded 
3-element airfoil with a short span of 5.7 cm (2.24 in) were performed with 218 million fine equivalent voxels. The 
current simulations were only intended to provide some guidance on how the model will perform aerodynamically 
and aeroacoustically while the detailed design could be influenced. However, calculations with finer grids are planned 
to make quantitative assessments of the noise from this model. Nonetheless, the HL-CRM will be constructed with 
removable slat edges around the pylon so that alternative geometries can be installed that may alleviate the noise source 
around the slat break. In the event that no reasonable modification can be found, the full-span slat case will be used to 
evaluate the slat noise reduction technology under development. 
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(a) Top view (b) Front view 


Figure 1. Two views of the HL-CRM geometry. 


Figure 2. PowerFLOW® grid distribution on a plane through the HL-CRM. 
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(a) Top (b) Bottom 


Figure 3. Identification of some of the variable resolution regions on the HL-CRM. 


(a) vOm Pylon (b) vOm Wing Tip 


(c) v2m Pylon (d) v2m Wing Tip 


Figure 4. PowerFLOW® grid distribution on planes through the HL-CRM. 
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(a) FUN3D Top (b) FUN3D Bottom 


(c) PowerFLOW® Top (d) PowerFLOW® Bottom 


Figure 5. Contours of the mean C,, on the HL-CRM. PowerFLOW® results from v2m mesh. 
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(a) FUN3D pylon (b) FUN3D outboard flap 


(c) PowerFLOW® pylon (d) PowerFLOW® outboard flap 


Figure 6. Contours of the mean C’, on parts of the HL-CRM. 


Figure 7. Planar cuts along the wing used for comparisons. 
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Figure 8. C> on planar cuts through the HL-CRM wing. 
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(a) Top view (b) Bottom view 
P 


Figure 9. Contours of the surface pressure fluctuations (p’,,,,, in Pa) on the HL-CRM. 


(a) B787 (b) A320 


Figure 10. Images of modern B787*’ and A320°° aircraft in the landing configuration. Photos obtained from Wikimedia Commons and 
reproduced under the Creative Commons Attribution-Share Alike 2.0 Generic license. 
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(a) HL-CRM (b) B787 (c) A320 


Figure 11. Comparison of the pylon / wing juncture region between the HL-CRM, B787,°” and A320.58 Photos obtained from Wikimedia 
Commons and reproduced under the Creative Commons Attribution-Share Alike 2.0 Generic license. 


(a) Baseline (b) Chine 


(c) Extended Slat (d) Full-Span Slat 


Figure 12. Views of contours of the surface pressure fluctuations (p/.,,., in Pa) on the HL-CRM. 


15 of 18 


American Institute of Aeronautics and Astronautics 


{/ 
NEN 
(c) Extended Slat (d) Full-Span Slat 


Figure 13. Close-up views of contours of the surface pressure fluctuations (p/,,,,, in Pa) on the nacelle/pylon of the HL-CRM. 


(b) Chine 


(d) Full-Span Slat 


Figure 14. Close-up views of contours of the surface pressure fluctuations (p/.,,., in Pa) on the pylon of the HL-CRM. 
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(a) Wing Tip (b) Top of Wing Tip 


(c) Outboard Flap Edge (d) Wing/Fuselage Junction 


Figure 15. Views of contours of the surface pressure fluctuations (p'.,.,, in Pa) on the HL-CRM. 


(a) Baseline (b) Chine 


(c) Extended Slat (d) Full-Span Slat 


Figure 16. Beamform maps for 2.5 kHz. 
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(a) Baseline 
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(c) Extended Slat 


Figure 17. Beamform maps for 5.0 kHz. 


(b) Chine 


(d) Full-Span Slat 


(a) Baseline 


(b) Chine 


(c) Extended Slat 


Figure 18. Beamform maps for 20.0 kHz. 
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(d) Full-Span Slat 


